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FOREWORD 

This  papfc-  was  prepared  by  the  Technical  Programs  Division, 
Department  of  Defense  Explosives  Safety  Board,  in  response  to  the 
need  for  techniques  for  rapidly  estimating  blast  and  fragment  hazards 
for  specific  an*. inition  storage  situations.  The  paper  draws  exten- 
sively upon  weoucnc  effects  data  compiled  by  the  Board  staff  and  its 
contractors,  in  Moated  by  the  literature  cited.  It:  is  hoped  that 
the  paper  will  be  of  wine  to  chose  concerned  with  analyzing  the 
potential  effects  of.  accidental  explosions  of  stored  ammunition  when 
normal  quantity-distance  standards  cannot  be  met. 


P.  F.  KLEIN 
Captain,  USN 
Chairman 
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PREFACE 

A computer  program  has  been  prepared  in  Fortran  language  for 
calculating  expected  numbers  of  casualties  by  blast  and  fragment 
injury  in  a population  of  unprotected  persons,  and  damage  to  un- 
strengthened residential  buildings  by  blast  from  explosions.  The 
program  generates  tables  of  injury  probability,  percent  casualties, 
and  percent  building  damage  as  functions  of  distance  from  block  stacks 
of  mass-detonating  bombs  and  projectiles  described  in  the  input  data. 

The  program  is  based  on  simple  representations  of  explosion  effects 
data  current  at  the  time  of  this  publication.  Although  injury  and 
damage  mechanisms  due  to  explosive  blast  are  considered  richer  well 
understood,  estimates  of  fragment  effects  are  based  on  extrapolation 
of  single-weapon  test  data  to  large  numbers  of  weapons  and  long  distances 
from  the  explosion  source.  The  code,  however,  can  be  readily  modified 
by  appropriate  substitution  of  improved  treatments  of  injury  mechanisms 
in  the  program,  or  better  characterization  of  weapon  properties  in 
input  data,  as  they  become  available. 

While  the  results  do  not  reduce  the  need  for  strict  compliance 
with  published  explosives  safety  standards,  they  may  be  of  assistance 
in  evaluating  the  protection  afforded  by  the  standards  and  in  assessing 
the  risks  connected  with  deviations  from  them  when  waiver  or  exemption 
requests  are  evaluated  or  granted. 
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COMPUTER  PROGRAM  FOR  PREDICTING  CASUALTIES 
AND  DAMAGE  FROM  ACCIDENTAL  EXPLOSIONS 


INTRODUCTION 


This  paper  describes  a computer  program  for  calculating  expected 
nurabe; s of  casualties  Trom  blast  and  fragments  in  a population  of 
u" protected  persons,  and  damage  .o  typical  un&trengthened  residential 
! ding;.  by  'last,  as  functions  of  distance  from  hypothetical  accidental 
c , <;-,ons  of  5ct  igular  stacks  of  mass-  Jetonating  bombs  and  projectiles. 

In  hat  folio  ■ < plosive  blast  and  fragmentation  effects  data  which 
form  the  basis*  ior  calculations  performed  by  the  code  are  summarized. 
Procedures  utilized  in  the  code  for  computing  probabilities  of  human 
injury  or  fatality  from  the  combined  effects  of  blast  and  fragments  are 
discussed.  Input  data  format,  program  opexation,  and  typical  results 
are  described.  A Fortran  listing  of  the  code  is  apoended,  together  with 
a sample  output  table  for  a particular  case. 

EXPLOSION  SOURCE  FUNCTIONS 

Blast  and  fragment  effects  are  calculated  from  the  equations  of 
carves  fitted  to  data  on  these  effects  as  functions  cf  distance  from 
the  source.  Blast  effects  are  referenced  to  a bare  TNT  hemisphere  on  a 
flat  ground  surface. 

Blast  Effects 

Blast  parameters  which  determine  human  injury  and  wood-frame  house 
damage  are  the  overpressure,  the  impulse  of  the  dynamic  (blast  wind) 
pressure  wave,  and  to  a minor  degree  explicitly,  the  wave  duration. 
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These  parameters  all  depend  on  distance  scaled  by  the  cube  root  of  the 
weight  of  the  equivalent  TNT  hemisphere.  Damage  to  certain  types  of 
buildings  such  as  large,  massive  structures,  which  may  be  a function 
of  the  impulse  of  the  overpressure  wave,  is  not  treated,  nor  are  close- 
in  (very  high  pressure)  effects  or  the  localized  effects  of  small  charges. 

In  the  computer  program  for  casualty  and  damage  prediction,  the 
effective  charge  weight  is  obtained  by  adjusting  the  total  net  explosive 
weight  by  factors  which  account  for  the  TNT  equivalence  of  the  explosive 
and  for  the  suppressing  effect  of  metal  case  material  and  magazine  earth 
cover,  if  either  or  both  are  present,  on  explosive  blast. ** 

The  case  effect  is  accounted  for  by  applying  a factor  F which  de- 
pends on  the  metal -to- charge  weight  ratio  M/C  in  the  following  way^: 

F = 0.6  - 0,4/U  * 2M/C) 

Ir  this  expression,  C is  taken  to  be  the  TNT  equivalent  of  the  filler 
explosive  weight  in  the  weapon. 

The  effect  of  the  earth  c-jver  of  a magazine  depends  on  the  cover- 
to-charge  ratio.  F ■ simplicity  the  single  numerical  value  of  0.8, 
derived  from  test  data  for  a standard  igloo  containing  ’00,000  lb  of 
explosive,  is  applied  to  the  explosive  weight  in  earth-covered  storage. 
Finally,  a factor  of  1.2  is  applied  to  account  for  the  enhancement  of 
blase  in  the  directions  perpendicular  to  the  faces  of  rectangular  stacks. 
This  value  has  been  found  to  be  a suitable  average  for  nearly  cubical 
stacks.* 

it 

Superscript  numerals  designate  appended  references 
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With  distance  R scaled  by  the  cube  root  of  the  weight  V of  the 
effective  TNT  hemisphere,  the  peak  overpressure  p is  calculated  from 
a ninth-order  polynomial  fit  of  its  logarithm  to  log(R/Wl/,3) , due  to 
Rosttsberg.3  Alternatively,  a similar  eighth-order  fit  given  by 
Kingery^  could  be  used. 

The  scaled  dynamic  overpressure  impulse  per  unit  area,  J/W*^3, 
is  represented  as  a one-term  power-law  function  of  scaled  distance 
adapted  by  Fugelso  et  al.*  from  Richmond  and  Fletcher.*’  Finally,  the 
scaled  positive  phase  duration,  t^/V^^,  is  also  expressed  as  a one- 
ten  power- law  function  of  scaled  distance  developed  by  Fugelso  et  al.* 

Fragment  Density 

Fragment  density  (number  per  unit  area)  is  calculated  as  a suitable 
multiple  of  that  from  the  sidespray  of  a single  weapon.  The  results  of 
Feinstein**  for  single  weapons,  obtained  as  a function  of  distance  by 
considering  the  effects  of  atmospheric  drag  and  gravity  on  fragment 
trajectories,  using  as  initial  data  the  velocities  and  weight  distribu- 
tions determined  in  weapon  effectiveness  tests,  serve  as  the  basis.  The 
single-weapon  density  profiles  so  determined,  for  all  fragments  and  for 
fragme  - of  energy  exceeding  58  ft-lb  at  impact,  are  represented  hy 
piccewise  linear  approximations  of  the  logarithm  of  fragment  density  as 
a function  of  distance.  The  approximation  for  heavy  fragments  consists 
of  two  straight  lines,  while  for  some  weapons  the  approximation  for  all 
fragments  requires  three  segments.  As  will  be  noted  later,  the  junction 
points  in  the  piecewise  linear  representations  are  identified  by  the  dis- 
tances at  which  specified  levels  of  fragment  density  occur  from  single  weapons. 
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In  applications  to  stacks  o£  weapons,  the  single 'weapon  density 
profiles  are  multiplied  by  the  number  of  units  on  the  side  or  end  of  the 
stack  considered  to  be  the  controlling  orientation.  No  distinction  is  made 
between  weapons  stored  aboveground  and  under  earth  cover.  Empirical  rules 
have,  however,  recently  been  suggested  which  might  be  incorporated  into 
the  code  to  improve  the  treatment  of  fragments  from  weapon  stacks. ^ 

INJURY  AND  DAMAGE  FUNCTIONS 

Personnel  injury  and  building  damage  are  computed  from  blast  pressure 
and  impulse  and  from  fragment  density,  all  as  functions  of  distance  from 
the  explosion.  Probabilities  of  injury  exceeding  each  of  two  levels  are 
considered:  severe  injury  or  fatality,  and  injury  of  any  severity,  however 

slight. 

Blast  Injury 

Human  injury  from  blast  may  be  one  of  three  types:  primary,  due  to 

the  overpressure  of  the  blast  wave  alone;  secondary,  from  impact  by  debris 
from  damaged  structures;  and  tertiary,  by  striking  the  ground  surface  or 
a stationary  object  after  acceleration  by  blast  forces. 

Damage  to  the  lungs  is  taken  to  be  indicative  of  primary  blast  injury 
of  a serious  level.  It  is  a function  of  blast  overpressure  alone,  except 
for  a weak  dependence  on  wave  duration.®  The  probability  of  severe  or 
fatal  injury  is  calculated  from  a cumulative  lognormal  distribution  with 
respect  to  the  pressure  (that  is,  from  a standard  error  function  of  the 
logarithm  of  pressure)  centered  at  £0.5  psi  for  waves  of  very  long  duration, 
with  1-  and  99-percent  levels  at  14.5  and  29.0  psi,  respectively. 
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Eardrum  damage  is  the  farthest- reaching  primary  effect  of  blast 
on  personnel.  It  is  a function  of  overpressure  only,**^  and  is  cal- 
culated from  a lognormal  <?  ^tribution  centered  at  6.3  psi,  with  10-  and 
90-percent  levels  at  3.2  and  12.2  psi. 

Secondary  effects  of  blast  on  personnel,  from  impact  by  building 
debris,  will  depend  strongly  on  the  extent  of  damage  to  structures.  As 
a result,  the  risk  to  occupants  can  at  best  only  be  inferred  qualitatively 
from  building  damage  calculations.  On  the  other  hand,  displacement 
(tertiary)  effects  can  be  estimated  for  persons  assumed  t.o  be  without 
protection  from  the  dynamic  pressure  impulse.*  Serious  head  injury  is 
assumed  to  occur  in  50  percent  of  cases  of  impact  at  18  ft/sec.  This 
corresponds  to  an  applied  impulse  or  75.2  psi-msec.  The  distribution 
is  normal,  with  1-  and  99-percent  levels  at  54.3  and  96.1  psi-msec.  A 
marginal  level  of  injury  by  this  same  mechanism  is  assumed  to  occur  in 
50  percent  of  exposures  2 ft/sec  (8.3  psi-msec),  with  a normal  dis- 
tribution having  1-  and  99-percent  levels  at  1.25  and  15.35  psi-msec, 
respectively.-* 

Fragment  Injury 

The  probability  P cf  a strike  by  one  or  more  fragments  is  calculated 


from  formulas  of  the  form 


P = 1 - exp(-qA) 


where  q is  the  fragment  density  and  A is  the  exposed  vulnerable  area  of 
the  individual.  For  calculations  of  severe  injury  probability,  the  density 
of  fragments  of  impact  energy  exceeding  58  ft-lb  is  used  in  the  formula, 
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with  area  A = 6.2  ft2,  corresponding  to  the  area  of  a person  facing  the 
explosion  and  taking  no  evasive  action.^  For  calculations  ox  all  fragment 


injury,  however  slight,  the  density  of  all  fragments  is  used  with  an  area 
of  8.0  ft2,  the  silhouette  area  of  a standing  man.9  The  former  value  of 


A is  considered  to  be  the  area  susceptible  to  serious  injury  by  heavy 


fragments,  while  the  latter  value  includes  body  extremities  as  well. 


Casualty  Probabilities 


Compound  probabilities  Pc  are  calculated  for  severe  or  fatal  injury. 


and  for  all  injury,  by  combining  the  probabilities  of  direct  blast  injury. 


displacement  injury,  and  fragment  injury  in  formulas  of  the  type 


P„  » 1 - Q-PtX1~P0)U-P„) 


where  Pj,  P2>  and  P^  are  the  probabilities  of  injury  by  the  independent 


mechanisms  operating.  The  result,  multiplied  by  100,  represents  the 


percentage  of  casualties  in  a population  exposed  at  a given  distance  from 


the  explosion. 


Building  Damage 


Damage  to  unstrengthened  single-family  dwellings  in  several  exposures 
to  test  explosions  has  been  assessed  quantitatively  by  Wilton.*0  He  found 


that  the  level  of  damage,  expressed  as  a percentage  of  building  value, 


depends  on  overpressure  alone.  The  results  are  well  approximated  by  a 


cumulative  normal  distribution  centered  at  3.8  psi,  with  15-  and  85-percent 


levels  at  2.1  and  5.5  psi. 
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COMPUTER  PROGRAM 

In  an  effort  to  provide,  in  easily  understood  tabular  fora,  overall 
blast  and  fragment  estimates  for  specific  ammunition  storage  situations, 
a computer  program  was  prepared  in  Fortran  language  incorporating  injury 
and  damage  functions  discussed  in  the  preceding  sections  of  this  paper. 

The  program  generates  tables  of  numbers  of  casualties  and  building  damage 
listed  as  functions  of  distance  from  100  to  10,000  ft  divided  into 
approximately  50  intervals. 

The  code  consists  of  a main  program,  two  function  subprograms,  and 
a block  data  subprogram.  Significant  steps  in  the  main  program,  identified 
by  line  numbers  in  the  appended  listing,  are  discussed  below. 

Line  2:  Storage  space  is  allocated  for  array  variables  containing 

descriptive  information  regarding  munition  stack  shape  (SS) 
and  the  format  of  numerical  output  (FMT). 

Lines  3-7:  Functions  are  defined  which  describe  the  logarithm  of 

fragment  density  as  piecewise  linear  with  respect  to  distance 
from  a single  munition. 

Lines  8-13:  Probability  units  are  defined  for  dynamic  overpressure 

impulse,  for  the  overpressure,  and  for  the  logarithm  of 
overpressure,  to  be  used  later  in  calculating  injury  and 
damage  levels  with  assumed  cumulative  normal  distributions. 

Lines  14-15:  A logarithmic  increment  of  distance  is  defined  fer 

later  use. 
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Line  17: 


Weapon  name,  gross  weight,  weight  of  exolosive  per  unit, 

TNT  equivalence  of  the  filler,  and  characteristic  fragment 
distances  from  a single  munition  are  read  in. 

Line  19:  Total  explosive  weight,  earth  factor  (1.0  or  0.8),  number 

of  units  on  the  end  of  the  stack,  stack  shape,  and  stack 
length  in  unit  diameters  are  read  in. 

Line  20:  If  stack  length  is  nonzero  in  input  data,  the  stack  is 

assumed  to  be  prismatic,  and  the  actual  total  filler  weight 
is  computed.  Otherwise  the  value  WE  supplied  in  input  data 
is  used. 

Lines  21-28:  Problem  definition  and  table  headings  are  written. 

Lines  29-31:  The  modified  Fano  case  factor^  and  the  effective  TNT 

weight  and  its  cube  root  are  calculated. 

Lines  32-33:  Distance  is  initialized  at  100  ft. 

Lines  34-37:  Scaled  distance,  overpressure,  duration  and  dynamic 

overpressure  impulse  are  calculate  i.  The  function  sub- 
program for  overpressure  is  based  on  a ninth-degree 
polynomial  fit.^ 

Lines  38-40:  Probabilities  of  serious  head  and  lung  injury  are 

calculated  from  cumulative  normal  and  lognormal  distributions. 
The  error  function  subprogram  is  based  on  a rational  approx- 
imation given  in  Abraroowitz  and  Stegun's  (1964)  NBS  Handbook 
of  Mathematical  Functions,^ 
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Line  41:  The  probability  of  blowdown  with  slight  injury  is  computed 

based  on  threshold  and  50-percent  levels  suggested  by 
Richmond. 

Line  42:  The  probability  of  ear  injury  (considered  not  a severe 

lesion)  is  calculated  from  a cumulative  normal  distribution 
with  respect  to  the  logarithm  of  overpressure. 

Line  43:  The  probability  of  primary  blast  injury  is  obtained  by 

compounding  ear  e n i lung  injury  probabilities. 

Lines  44-56:  The  density  (Q)  of  fragments  of  impact  energy  greater 
than  58  ft-lb,  and  the  density  (QA)  of  all  fragments,  are 
calculated  unless  the  fragment  distance  parameter  X3  has 
not  been  supplied  as  input.  Q is  defined  as  bilinear  in 
semilogarithmic  coordinates,  and  QA  as  either  bilinear  or 
trilinear,  coinciding  with  Q at  distances  greater  than  X3. 

The  piecewise  linear  approximations  for  the  155mm  projectile 
with  base  down,  shown  in  Figure  1,  are  typical.  At  the 
distances  XO,  XI,...,  X4  supplied  in  input  data,  the  fragment 
density  from  a single  weapon  has  values  which  are  specified 
to  be  multiples  of  d/60Q)ft"^  by  factors  of  10*^  in  the 
arithmetic  statement  function  definitions,  lines  3-7. 

Lines  58-59:  The  probabilities  of  impact  by  heavy  fragments  end  by 

all  fragments  on  the  presented  area  of  an  individual  facing 
the  explosion  are  calculated  from  Poisson  statistics.  The 


9 


former  represents  severe  injury  or  death,  while  the  latter 
is  an  upper  bound  intended  to  represent  all  fragment  injury, 
however  slight. 

Lines  60-61:  Serious  injuries  per  100  population,  and  all  injuries, 

are  calculated  by  compounding  the  appropriate  probabilities. 

Line  62:  Damage  in  percent  of  building  value  is  calculated  from  a 

cumulative  normal  distribution  with  respect  to  overpressure 
centered  at  3*8  psi.  This  combines  structural  and  super- 
ficial damage,  and  is  an  excellent  fit  of  data  compiled  by 
Wilton.10 

Line  63:  Results  are  written  in  formats  specified  at  program  execution 

time  (line  16).  The  formats  UHO,  F8.0,  F9.2,2(3F10.2, 

F10.0) ,F10,0)  and  <1HO,F8,0,4F10.3,F10.1 ,3F10.3,2F10.1) 
have  been  used  for  2-  and  3-place  precision  in  output. 

Lines  64-69:  The  running  distance  variable  X is  increased  by  a 

constant  logarithmic  increment,  and  R is  the  next  lower 
integer  multiple  of  10  ft,  not  exceeding  the  table  range  of 
10,000  ft.  A coarser  table,  with  values  of  R being  multiples 
of  100  ft,  is  obtained  if  one  replaces  line  67  with  the  statement 
R = 100*1FIX(X/100. ) . 

A sample  table  of  output  data  at  2-place  precision  is  appended  with 
the  program  listing.  Two  data  cards  are  required  to  generate  each  such 
table.  In  addition,  the  first  data  card  of  the  deck  must  specify  output 
format  for  an  entire  job. 


The  tabulated  results  may  be  used  directly  if  desired,  or  they 
may  be  plotted  and  used  graphically.  Curves  of  injury  and  damage  for 
approximately  250,000  lb  of  tritonal  in  750-lb  bombs  are  shown  in 
Figure  2.  The  results  of  a number  of  cases  may  be  cross-plotted  to 
determine  the  relations  between,  say,  explosive  quantity  and  distance. 
An  example  is  given  in  Figure  3,  where  such  curves  are  shown  for  two 
types  of  weapons  at  the  10  percent  level  of  serious  injury  or  fatality. 


■ 
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PROGRAM  CASDAMC INPUT » OUTPUT »TApF60a INPUT, TApE6 l *UUTPUT) 
DIMENSION  3S{2) pFMTCB) 

0I2(X)  * 10,**(2,*(X»X2)/(3,*(Xl«X2))«8,/3()/6, 

U2i(X)  • 10»**(  <XaX2)/(3,*<X2*X3))aB,/3,)/6, 

5 Q03CX)  « 10,**C5,*(X»X3)/C3,*Cx0*X3>)«*3.)/6, 

80PCX)  c io,**<4,*<x*X4)/c3,*cxe*X4))»8,/5,)/6, 

Q*i<x)  • io,**c  Cx*X4)/< 3t*(xa-X3)}*e,/3.)/6f 

3QRT2  ■ SQRTC2.) 

SOJ2  • SQRT2  *41, 8/4,653 

10  3DJ2T  • S0RT2  *7,06/2,327 

SOI.P2  » 80RT2  *AL0G(29, 0/14, 5)/4, 653 

30LP2E  a 3QRT2  *AL06(12, 2/6, 31/1, 282 

SOP2  • SQRT2  *3,4/2,073 
3PC  * 25, 

15  UU1  * 10,**(1 ,/SPC)+l ,E«7 

RfcADC60,100)  FMT 

1 READC60, 101)  MPK,MGR,rtEU,T£Q,Xl ,X2;X3,X0,X4 
IF CE()F (60),NF,0,)  STOP 

READC69, 102)  *E,EF,B3,33,3* 

20  1F(3*,GT,0,)  a hS*SW*neU 

iw  ■ Nt 
18  ■ 83 

IE  a S,«4,99*EF 

wRITE(61 ,201 ) *PN,IW,TEQ,IB,IE,SS 
25  HRITE(6i,202) 

HR1TE(61,203) 

MRITllol ,204) 

NRITfc(6l,205) 

FF  ■ 0,6*0,4/C!,*2,*{NSR»NEU)/C1EO*WEU)) 

30  at/  « l,2*NE*TEi)*FF*fcF 

CRT*  a wEF**{!./3,) 

X a 100, 

Rax 

2 L a R/CRTN 

35  PR  a PRF(Z) 

TO  a 1 , 188*CRTm*Z**0, 3473 

XJ  a 266,449*CRTP/Z**2,320 1 

PHO  a (l,*ERF<CXJ*75,2)/3DJ2))/2, 

P«TD  a 20,5*(1 ,*6, I6/TD**1 ,064) 

40  PUG  a U,tFRF(AL0G(PR/PMTD)/SDUP2))/2. 

PIC  a (l,*ERFt(XJ-8,3)/SDJ2T))/2, 

PtR  * Cl ,*ERFCALOG (PR/6,3) /SDUP2E))/2, 

P8I  ■ l,*U,*>PER)*Ci,«PLC) 
tf  ■ 0, 

45  OA  a 0, 

IF <Xi,UE#0, ) CO  TO  3 
IFCR.LT.X2)  Q a BS*Q12(R) 

IF (R,6E ,X2)  G a 8S*Q23(R) 

UA  ■ Q 

50  IFCR,GE,X3)  GO  TO  3 

IF (X4,GT,0,)  GO  TO  5 
UA  a BS*Q03(R) 

GO  TO  3 
6 CONTINUE 

55  IFCR.LT.X4)  QA  a B3*O04CR) 

IFCR.CE.X4)  QA  a BS*943(R) 

3 CONTINUE 
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70 
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60 


65 
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PFK  • 1,«>EXP(*6,2*Q) 

P*P  ■ 1,«£XP<-8,0*OA) 

PCTC  « 100,*d,«d,-PHD)*d.»PLG)*d,*PFR)) 

PCTAC  • lOO.^d.od.-PIGJtktij-PBIJod.tiPAP)) 

PCTO  ■ 100,*d,  + ERP((P«»3,6)/80P2))/2, 

«RIT£(61,FMT)  H,Z«PHD,PLG#PFR»PCTC,PI0,P8I#PAP, PCTAC, PCTO 
Hu  # R 
4 CONTINUE 
X « Dl.X  1 *X 
R • 10*IF!XJX/10,> 

IF(R.LE.RO)  GO  TO  4 
IF{R.GT,1,E4)  GO  TO  I 
GO  TU  2 

100  FORMAT  ( 6 A 1 0 ) 

101  FORMAT  (A6,6F10,0,2P6,OJ 

102  FORMAT  <3F10,0,2A8,F10,0> 

201  FORMAT  (1M1,A8»U0  #9M  LB  ME  ATF5a2,8H  TNT  £0,14  »24H  SIDE  UNITS 

1»  EARTH  C0VERI2#14H,  STACK  SHAPE  2A6) 

ItS^OFlei j72pfR?rif?,*RC1B*B1LlTIE8  ^ *®*» ImPERCFNT#6X»  16HPR0BABTLITI 

2°34^nDiTr,<1H  *25X#22HSEVERE  OR  FATAL  lNJURY#5X#7M3ERlfVJSa6X»22MANY  I 
INJURY  OR  FATALITY, 7X , 3HANY,SX , 7HPERCENT > 

204.^Hp*T  ‘^'“ROi^.RX^HSCAUEOrSXfRMHEAO^XfKHLUNG.fcX^MMEAVY, 
2L01NG)  XNjURY'3x'6HGROU''IO'!>*'5HBLAST»6X,3HALL*6X,6HlNJURY,2Xa6H6UI 

2°5irM^rJ  lJH.ifX'^^T>'5x»,MDIST'5x*fcMlN'»URY*3’<»flHCOLLAPSr»2X,4MFRA 

2R  DtATMf2x!lHDA25J^'2X'6HlMPACT'aX,*MINJURV,5X',MFRAGMtNTS'2X#8Hn 
ENU 
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FUNCTION  PRF(Z) 
COMMUN/OV/OP(IO) 

Zt  a ALOCCZ) 

SUN  a OP(tO) 

DU  t I i l#4 

IN  B m\ 

i SUN  a 9tJM*ZL+OP(INMO) 
P«F  * EXP(8UH) 

RETURN 

END 


FUNCTION  ERF (X) 
CON*UN/£R/A(5) 

AX  b X 

IFCX.LT.O.)  AX  a .X 
SUN  a A (5) 

DO  1 I a 1,4 
IN  a ml 

l SUN  ■ 8UM*AX*A(IN45) 
ERF  « l.«la/SUN**4 
IF(XjLT,0#)  ERF  • *ERF 
RETURN 
END 


SLOCK  DATA 
CONNON/OV/OP(»0) 

CONNQN/ER/ATS) 

data  OP/,705fcBlEi,.,l6S37ZAEl,*.25l6«ei,-,itJ771«.,T818AOSE-l. 

I ,S0J5lR8E-J,.t275*R7E-l,,5557R*ef*2,-,5i08AI«E-S>,17R55*5fc»«/ 
DATA  A/l (# ,278395* ,230389*  » OS 9972# ,078108/ 

END 
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COMPUTER  PROGRAM  FOR  PREDICTING  CASUALTIES 
AND  DAMAGE  FROM  ACCIDENTAL  EXPLOSIONS 


Tables  of  Results 


Numerical  results  tabulated  on  succeeding  pages  were  calculated 
by  means  of  a computer  program  described  in  Department  of  Defense 
Explosives  Safety  Board  Technical  Paper  11,  of  which  this  is  a 
supplement.  The  program  generates  tables  of  injury  probability, 
percent  casualties,  and  percent  damage  to  residential  buildings  as 
functions  of  distance  from  block  stacks  of  mass-detonating  bombs  and 
projectiles  described  in  the  input  data. 

Results  are  presented  for  aboveground  stacks  of  the  following 
four  weapons:  500- lb  bomb  Mark  82,  H-6  loaded;  750- lb  bomb  M117, 

tritonal  loaded;  155-ram  projectile  M107,  TNT  loaded;  and  175-ram 
projectile  M437A2,  Composition  B loaded.  The  stacks  considered  are 
assumed  to  be  built  up  of  weapons  strapped  to  their  normal  shipping 
pallets,  and  range  from  approximately  5,000  to  500,000  lb  in  net 
explosive  content.  In  these  configurations,  the  bombs  are  stored  with 
axes  horizontal,  while  the  projectiles  are  stored  with  bases  down. 

The  scaled  distance  listed  in  the  second  column  of  each  table  is  in 
units  of  feet  divided  by  the  cube  rcot  of  the  equivalent  weight,  in 
pounds,  of  uncased  TNT. 
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hark  02  5000  Ht  AT  1,2b  TNT  EO»  0 SIDE  UNITS,  EARTH  COVER  0,  STACK  SHARE  2 


M117  500000  lb  HE  AT  1«13  TNT  £0/  SB  SIDE  UNITS#  EARTH  COVER  0,  STACK  SHAPE  12  X SB  X S 

probabilities  OF  percent  probabilities  op  percent 
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probabilities  of  percent  probabilities  of  percent 

SEVERE  OR  FATAL  INJURY  SERIOUS  ANY  INJURY  OR  FATALITY  ANY  PERCENT 
DIST  SCALED  HEAD  LUNG  HEAVY  INJURY  GROUNO  BLAST  ALL  INJURY  BUILDING 
(FT)  OIBT  INJURY  COLLAPSE  FRAGMENTS  OR  DEATH  IMPACT  INJURY  FRAGMENTS  OR  OeATH  DAMAGE 


Mil?  SOJOO  LB  Mfc  A r 1,15  TNT  F.8,  12  $10£  UMTS,  EARTH  COVER  0,  STACK  SHAPE  4 X 12  X 3 

PROBABILITIES  OF  PERCENT  PROBABILITIES  OF  PERCENT 

SEVERE  OR  FATAL  INJURY  SERIOUS  ANY  INJURY  OR-  FATALITY  ANY  PERCENT 

OlST  SCALEO  Hfc.Au  LUNG  HEAVY  INJURY  GROUND  BLAST  ALL  INJURY  BUILDING 


PROBABILITIES  OF  PERCENT  PROBABILITIES  OF  PERCENT 

SEVERE  OR  FATAL  INJURY  SERIOUS  ANY  INJURY  OR  FATALITY  ANY  PERCENT 

OIST  SCALED  HEAD  LUNG  HEAVY  tNJORY  GROUND  BLAST  ALL  INJURY  UUILOINO 

(FT)  OIST  INJURY  COLLAPSE  FRACHENT9  OR  DEATH  IMPACT  INJURY  FRAGMENTS  OR  OEATH  DAMAGE 


t, 


1*5  HN  250000  u»  H£  AT  1,00  TNT  £0*  164  8J0E  UNITS,  EARTH  COVER  0,  STACK  SHARE  68  X 92  X 2 

PROBABILITIES  of  percent  probabilities  op  percent 

SEVERE  OR  FATAL  INJURY  SERIOUS  ANT  INJURY  OR  FATALITY  ANY  PERCENT 

OIST  SCALED  HEA0  LUNG  HEAVY  INJURY  GROUND  BLAST  ALL  INJURY  BUILDING 

(FT)  OIST  INjUNy  COLLAPSE  FRAGMENTS  OR  DEATH  IMPACT  INJURY  FRAGMENTS  OR  DEATH  DAMAGE 
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t*5  MM  5000  L»  HE  AT  1.00  TNT  EG,  29  8XOE  UNITS,  EARTH  COVER  0,  STACK  SHARE  tU  tU  2 

PRO“ieZUTIES  OP  percent  probabilities  of  percent 

SEVERE  OR  FATAL  INJURY  SERIOUS  ANY  INJURY  OR  FATALITY  ANY  PERCENT 

OI3T  SCALED  HEAD  LUNG  HEAVY  INJURY  GROUND  BLAST  ALL  INJURY  BUXLOING 

(FT)  GIST  INJURY  COLLAPSE  FRAGMENTS  OR  DEATH  IMPACT  INJURY  FRAGMENTS  OR  DEATH  DAMAGE 
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